Microstructure of a liquid complex (dusty) plasma under shear 
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The microstructure of a strongly coupled liquid undergoing a shear flow was studied experimen- 
tally. The liquid was a shear melted two-dimensional plasma crystal, i.e., a single-layer suspension 
of micrometer-size particles in a rf discharge plasma. Trajectories of particles were measured using 
video microscopy. The resulting microstructure was anisotropic, with compressional and extensional 
axes at around ±45° to the flow direction. Corresponding ellipticity of the pair correlation function 
g(r) or static structure factor S(k) gives the (normalized) shear rate of the flow. 

PACS numbers: 52.27.Lw, 52.27.Gr, 82.70.Dd 



Shear flows are ubiquitous components of almost any 
flow in liquids. They are important in many indus- 
trial applications and pose interesting questions in fluid 
mechanics. One question, which goes beyond the con- 
ventional hydrodynamic approach (based on the Navier- 
Stokes equation) is a relation between the strain rate in 
a flowing liquid and its microstructure. It has long been 
recognized that the microstructure of various sheared liq- 
uids is distorted in a way similar to an elastically de- 
formed solid (glass) In fact, the very existence of 
shear viscosity in a liquid is closely related to this dis- 
tortion. Experimentally, this distortion can be indirectly 
seen in the static structure factor of a sheared liquid de- 
duced from scattering measurements. Direct observation 
of the microstructure distortion in real space is impracti- 
cal for regular liquids and is only possible in simulations 
or model systems where the motion of individual particles 
- proxy "atoms" - can be resolved. Suitable model sys- 
tems are colloidal dispersions and complex (dusty) 
plasmas @, Q- 

Complex plasma is a suspension of fine solid particles 
in a weakly ionized gas d, 0j- Micron-size particles ac- 
quire high electric charge (usually negative) and may self- 
organize in ordered structures. The particle dynamics is 
governed by interparticle interactions as well as by in- 
teractions with the surrounding plasma. The latter can 
lead to various instabilities resulting in anomalous "heat- 
ing." However, for single-layer suspensions these instabil- 
ities are well understood and can be easily suppressed by 
proper choice of experimental parameters Q . In this case 
the particle kinetic temperature is a well-defined value 
determined by the temperature of ambient gas. Depend- 
ing on the particle coupling strength (i.e, the ratio of 
the pair interaction energy to temperature), a complex 
plasma can be in a liquid or solid state. The motion of 
individual particles is fully resolved in real time, allowing 
direct observation of their dynamics at the "atomistic" 
level. In two-dimensional (2D) complex plasmas, parti- 
cles are arranged in a single layer and are therefore easy 
to observe using video microscopy. 

Complex plasmas belong to a broad class of soft con- 
densed matter: Typical shear modulus of a crystalline 



FIG. 1: Shear flow in a planar Couette configuration in a 
strongly coupled 2D complex plasma is sustained by a pair of 
counterpropagating laser beams (indicated by arrows) . Parti- 
cle trajectories during 0.83 s of a steady-state flow are shown. 
The dashed-line box defines the region of interest where the 
analysis of liquid microstructure was performed. The inset 
shows the transverse profile of particle velocity v v (x). The 
manipulation laser power was Pi aS er = 1-25 W. 



2D sample is around 1CP 13 N/mm so that they can 
be easily manipulated, e.g.. by applying radiation pres- 
sure of a laser beam ljl Hi ■ This allows one to study a 
broad range of phenomena from nucleation and motion 
of single dislocations 
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to shear flows 

the shear viscosity of liquid complex plasma was calcu- 
lated using a fit of measured particle velocity profiles to 
a Navier-Stokes model. However, any changes occurring 
in the microstructure of sheared liquid complex plasma 
have not been studied so far. 

In this paper, we experimentally study the microstruc- 
ture of a sheared 2D complex plasma at various levels of 
strain rate. We find that the microstructure is distorted, 
with compressional and extensional axes at around ±45° 
to the flow direction. The magnitude of distortion is de- 
fined by the (normalized) shear flow's strain rate. 

Our experimental setup was a modified GEC (Gaseous 
Electronics Conference) rf reference cell [9j . Plasma was 
produced using a capacitively coupled rf discharge in ar- 
gon at 0.66 Pa. A single layer of dust particles was sus- 
pended in the plasma sheath of the lower rf electrode. 
The microspheres made of melamine formaldehyde had a 
diameter of 9.19 ± 0.09 fim, a mass to = 6.15 x 10~ 13 kg, 
and acquired an electric charge of Q = —17 000 ± 1700e. 
The suspension included around 8000 particles and had 
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FIG. 2: Angle-averaged (a) pair correlation function g(r) and 
(b) static structure factor S(k) of the complex plasma under 
different levels of shear stress as controlled by the manipu- 
lation laser power Pi asor (indicated in the inset). S(k) was 
calculated (in this figure only) in a central square of twice the 
width of the rectangular region of interest, see Fig. [T] 



a diameter of ~ 60 mm, the mean interparticle distance 
in the center was A = 0.55 mm [measured from the first 
peak of the pair correlation function g{r)]. The neutral 
gas damping rate was v = 0.77 s _1 . 

The like-charged particles repelled each other via 
a screened Coulomb, or Yukawa potential: U(r) — 
Q(47reor)~ 1 exp(— r/Xp), where Xp is the screening 



length 13]. The screening parameter n — A /Xd was 



around unity in our experiment |14| . The particles were 
held together by electrostatic fields naturally present in 
plasma. Since the vertical confinement of particles is 
much stronger than their horizontal confinement Q , they 
remained in a single layer at all times. We verified this 
by a side-view camera. 

Under our experimental conditions, the particle sus- 
pension self-organized in a highly ordered triangular lat- 
tice. We used the method of Ref. 12[ to create a shear 
flow in the 2D complex plasma as shown in Fig. [T] Two 
oppositely directed laser beams were focused down to a 
fraction of the interparticle spacing and they were rapidly 
(~ 300 Hz) scanned to draw rectangular stripes on the 
suspension. The particles reacted to the averaged ra- 
diation pressure. Shear stress was created in the gap 
between the laser-illuminated stripes, its magnitude was 
controlled by varying the output laser power. 

To facilitate the emergence of shear flow and subse- 
quent analysis, one of the closely packed rows of the tri- 
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FIG. 3: (a) Interparticle distance A and (b) bond density as 
functions of azimuthal angle <j> in liquid complex plasma at 
different levels of shear stress (controlled by Pi aser , see inset). 
The flow direction is indicated by the arrows in (a). The 
liquid's microstructure is anisotropic, with the compressional 
and extensional axes at around ±45° to the flow direction. 
The 60° periodicity of the bond density in (b) indicates that 
the liquid complex plasma retains to some extent the original 
hexagonal structure of the plasma crystal (see also Fig. 
The curves for unmelted crystal (P ase r < 0.9 W) are discon- 
tinuous and are not shown here. 



angular lattice was always oriented along the laser beams, 
before the laser was switched on. Unless otherwise stated, 
data analysis was performed at the stage of steady-state 
shear flow (after a waiting time of 3.3 s) in a region of 
interest indicated by a dashed-line rectangle in Fig. [T] 
All values reported were averaged over another 3.3 s. 

The transverse profiles of particle velocity in the 
steady-state shear flow, v y (x), were close to linear, see 
inset in Fig. [1] Therefore, we assumed that the shear 
strain rate 7 was uniform over the flow and calculated 
it as the slope of the linear fit of v y (x). In our system, 
the profiles of v y (x) should in fact be exponential (due to 
the presence of neutral gas drag force on particles), with 
the spatial scale defined by the (square root of) ratio of 
kinematic viscosity to the gas drag rate, see Ref. [l2| . In 
the present experiment, however, the width of the flow 
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was smaller (roughly two times) than in Ref. [12j, so 
the neutral gas friction played practically no role here. 
The transverse profiles of particle kinetic temperature T 
were relatively flat within the region of interest shown in 
Fig. [TJ with a variation of 8% — 30% depending on the 
manipulation laser power. In these conditions, the dif- 
fusive term in the heat balance for particles was at least 
one order of magnitude smaller than the friction term, 
so that the heat balance reduced to the following equa- 
tion: vnks{T — To) = rjj 2 /A, where n, r\ and To are, 
respectively, areal number density, shear viscosity, and 
background temperature of the particle suspension flBj ] . 
Experimentally measured temperature scaled roughly as 
T ~ Aj, where A = 5.7 x 10 4 Ks. This is consistent 
with the rj oc T -1 scaling in the above equation, in the 
relevant temperature range fl6j ] . 

We begin with standard structure analysis, i.e., calcu- 
lating the pair correlation function g(r) and static struc- 
ture factor S(k) of complex plasma (averaged over az- 
imuthal angle (/)), shown in Figs. EJa) and (b), respec- 
tively. Both indicate that the complex plasma is in a 
crystalline state for the laser power Pi ascr < 0.9 W and 
becomes liquid for Pi aser > 0.9 W, as evidenced by a sud- 
den decrease in the height of the first peak of g(r) and 
S(k) and disappearance of the splitting in their second 
peak. Although further subtler changes in both func- 
tions occur with increasing levels of shear strain rate in 
the liquid state, these are less prominent and difficult to 
decipher. 

To get a greater insight on the liquid microstructure, 
we propose to calculate the mean intcrparticle spacing A 
as a function of the azimuthal angle <f>. This is equivalent 
to evaluating the first peak in the angle-resolved pair cor- 
relation function <?(r), yet provides a very local measure 
and is also easier to calculate. Experimentally observed 
A((/>) is shown in Fig. |3](a) . A complementary measure 
is the angular bond density, i.e., the probability to find a 
near neighbor at a certain angle <f> from a given particle, 
see Fig. |3]». 

The microstructure of sheared liquid is clearly 
anisotropic: The interparticle spacing is larger in a cer- 
tain direction (extensional axis at 4> ~ —45°) and smaller 
in the orthogonal direction (compressional axis at <f> ~ 
45°). This constitutes a departure from isotropic distri- 
bution in a quiescent liquid. 

The distorted pair correlation function is given by 



g(r) = g (r)[l + T^{xy/r 2 )f{r% 



(1) 



where <?o( r ) is for liquid in equilibrium, 7 is shear rate, r 
relaxation time, and the function f(r) is determined by 
the particular form of the pair interaction potential for 
particles. The term proportional to xy/r 2 introduces a 
sin(20) component in the angular dependence of g(r). 
Similarly, for the distorted static structure factor one 
has S(k) = Sb (AM 1 + Tj{k x k y /k 2 )T(k)}. We note that 
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FIG. 4: Fits with Eq. © of experimental A(<f>) for (a) Pi aser 
1.05 W and (b) 1.2 W. 



these expansions for "strained" g(r) and S(k) are valid 
for small rj. 

In our experiment, the biggest contribution to the os- 
cillatory part of A(</>) is given by sfn(20), with a coeffi- 
cient A2 ~ 0.035. This term, however, does not explain 
all features of the observed A(</>) and the inclusion of 
higher harmonic sin(40) is necessary with A\ ~ 0.015. 
The fits of experimentally measured interparticle dis- 
tance with the following ansatz: 

A{(j)) = A [l + A 2 sin(20 + S 2 ) + Atsin(40 + S 4 )} (2) 

are shown in Fig. [4j 

The observed anisotropy of sheared liquid's microstruc- 
ture has clear physical meaning. The flow's strain rate 
(normalized) is given by the degree of anisotropy, i.e., 
rj ~ 2^2- The deviation of the directions of compres- 
sional and extensional axes from ±45° gives the extent 
of the non-Newtonian behavior of liquid The sig- 
nificance of A^ is less clear, it may reflect some kind of 
nonlinearity. 

Another way of evaluating the microstructure of a 
sheared liquid is to measure its static structure factor 
S(k) with angular resolution. For regular liquids, this 
is usually deduced from scattering experiments; here, we 
calculate S(k) directly from particle positions 17]. It is 
shown in Fig. [5] for different levels of applied shear stress. 
At very low stress, S(k) retains the structure of sharp 
peaks characteristic of an undisturbed triangular lattice, 
as shown in Fig. 0Ja). The peaks become blurred at 
higher stress as the lattice deforms and melts and even- 
tually they coalesce into rings in fully developed shear 
flow, Figs. EKc),(d). The rings are distorted from the cir- 
cular shape pertinent to quiescent liquids: They acquire 
certain ellipticity. This is particularly clearly seen in the 
first Debye-Scherrer ring. The ratio of its major and mi- 
nor axes is given by 1+e, where e = rj [HQ. This allows 
one to derive T7 from shear-distorted S(k). 

Figure [6] shows experimentally measured e and ^2,4 as 
functions of shear rate 7. There is a trend for e and A2 
to increase and for A4 to diminish for larger 7. We note 
that e is larger than 2^2 by a factor of 1.2 — 1.6. The 
uncertainty of e is relatively large, since e was directly 
measured from figures like Fig. [5] (not from a fit). 
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FIG. 5: Static structure factor S(k) of the 2D complex plasma 
under different levels of shear stress (controlled by Piascr, as 
indicated). The arrows in (d) indicate schematically the di- 
rection of the flow. Color coding from dark blue to red corre- 
sponds to the range of — 10 in (a) and — 5 in (b),(c), and 
(d). 



The resulting values of ry depend only weakly on 7 
and saturate, for larger 7, at about 0.11 and 0.07 for the 
5(k) and A(<f>) data, respectively. The relaxation time r 
is therefore not constant but scales roughly as r tx y^ 1 . 
(The inset in Fig. [6] shows that fits with r = const do not 
work well.) Recalling that the particle kinetic tempera- 
ture T ~ Ay, we then find that r ~ 0.1.4T -1 . In the 
range of T = 2 x 10 4 — 6 x 10 4 K achieved in the present 
experiment, r = 0.1 — 0.3 s. The temperature scaling 
of relaxation time hints at the particle diffusion as the 
process that defines it. We estimate the diffusion time in 
our system as T^ig — A 2 /D, where D is the self-diffusion 
coefficient. Because of the lack of systematic data on 
diffusion in two-dimensional complex plasmas, we use 
the empirical approximation for D in three-dimensional 
equilibrium Yukawa systems proposed in Ref. (lij : D ~ 
0.01cjEA 2 (T/T m — 1), where lue is the Einstein frequency 
[l9| and T m is the (equilibrium) melting temperature that 
should be regarded as a dimensional parameter in our 
case. Since the measurements of microstructure were per- 
formed well in the liquid state, we assume that T 3> T m 
and obtain wet^s ~ 10 2 T m /T, which is consistent with 
the experimentally observed temperature scaling of re- 
laxation time. Quantitative agreement with experiment 
is achieved at T m ~ 10~ 3 Aje = 3.2 x 10 3 K; this is much 
lower than the experimental range of temperature. 

To summarize, we studied for the first time the mi- 
crostructure of a 2D liquid complex plasma subjected to 
shear flow. To evaluate the microstructure, we used two 
measures: the mean interparticle distance as a function 
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FIG. 6: Ellipticity e of the liquid static structure factor 
S(k) (triangles) and relative amplitudes A2 (circles) and A± 
(squares) in Eq. ((2j), for different levels of shear strain rate. 
The inset shows linear fits (forced to 0) of e and A2 for 
7 < 0.8 s" 1 . 



of azimuthal angle A(</>) and the static structure factor 
5(k). Both methods reveal anisotropic microstructure, 
with compressional and extensional axes at around ±45° 
to the flow direction. The ellipticity of A(</>) or S(k) 
yields the (normalized) shear rate of the flow. The two 
measures give comparable results. 
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